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Abstract 



We calculate one-loop correction to the axial Ward-Takahashi identity given 
by Furman and Shamir in domain-wall QCD. It is shown perturbatively that 
the renormalized axial Ward-Takahashi identity is satisfied without fine tuning 
and the "conserved" axial current receives no renormalization, giving Za = 1- 
This fact will simplify the calculation of the pion decay constant in numerical 
simulations since the decay constant defined by this current needs no lattice 
renormalization factor. 
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Typeset using REVT^ 



I. INTRODUCTION 



The lack of chirally invariant fermion formulations is one of the most uncomfortable 
points theoretically and practically in lattice QCD. For example, in the Wilson fermion 
formulation, which is popularly used in numerical simulations, the chiral limit can be realized 
only by the fine tuning of bare mass parameter, which compensates the additive quantum 
correction to the quark mass. 

Recently the domain- wall fermion formulation which was originally proposed for 

lattice chiral gauge theories 0, has been employed in lattice QCD simulations [§| and has 
shown its superiority over other formulations: there seems no need of the fine tuning to 
realize the chiral limit while there is no restriction to the number of flavors. In particular 
massless mode which presents at the tree level seems stable against the quantum correction. 
This property of the domain-wall fermion seems to suggest the existence of a kind of axial 
symmetry. In Ref. Furman and Shamir have defined the axial transformation and the cur- 
rent in the domain-wall fermion. They have shown that although their transformation is not 
a symmetry of the action, the axial Ward-Takahashi identity is satisfied in Green's functions 
between the current operators and physical quark fields made of boundary fermions. 

Although it has already argued in the same reference that the axial Ward-Takahashi 
identity seems to hold nonperturbatively, it is still interesting to see a definite form of one 
loop correction to Green's functions in the identity since renormalized form of the Green's 
functions is nontrivial. In this paper we calculate the one loop correction to the Green's 
functions with the axial current operators and the quark fields. It is shown that all the 
divergence are renormalized into the same Green's functions without any mixing between 
these physical operators and unphysical heavy fermion operators. The renormalized axial 
Ward-Takahashi identity is satisfied without any fine tuning and the renormalization factor 
of the axial current becomes unity, Zj^ = 1 . 

This paper is organized as follows. In Sec. II we introduce the action of domain-wall QCD. 
The axial transformation and its Ward-Takahashi identity is given according to Ref. 0. In 
Sec. Ill we calculate one loop correction to the Ward-Takahashi identity. The Feynman rules 
relevant to our calculations are defined in our previous paper and are briefly given in the 
appendix. Our conclusion is given in Sec. IV. 

In this paper we set the lattice spacing a = 1 and take SU {Nc) gauge group with the 
gauge coupling g. 

II. AXIAL CURRENT AND WARD-TAKAHASHI IDENTITY 

The domain- wall fermion is a 4-1-1 dimensional Wilson fermion with a "mass term" which 
depends on the coordinate in the extra dimension. For the explicit form of the "mass term" 
we adopt the Shamir's one and the corresponding action becomes 

^ [l — — 

+ i (^{n),{l + 75)z^(n),+i + - 75)^(n),_i) + (M - 1 + 4r)V^(n),V^(n), 
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+ mJ2['ilj{n)NP+^in), + V^(n)iP_V^(n);v) , (II.l) 

n 

where n is a 4 dimensional space-time coordinate and s is an extra fifth dimensional or flavor 
index, the Dirac "mass" M is a parameter of the theory which we set < M < 2 to realize 
the massless fermion at tree level, m is a physical quark mass, and the Wilson parameter 
is set to r = — 1. It is important to notice that we have boundaries for the flavor space; 
1 < s < N. P± is a projection operator P± = (1 ± 75)/2. 

The remarkable property of the domain-wall fermion is that there exists a massless 
fermion mode in the N ^ oo limit at small momentum at m = 0. This massless fermion 
stays near the boundaries of the flavor space with the left and the right mode on the opposite 
side. At tree level the massless mode Xo is given explicitly in zero momentum limit as 



1 - w^o (P+wl^i^s + P~w^~'i^s] , (II.2) 



Xo = 

where wq = 1 — M. Although the zero mode is stable against quantum correction 0, the 
damping factor wq is renormalized due to the additive correction to the Dirac mass M. Thus 
in numerical simulations it is more convenient to use the interpolating "physical" quark field 
defined by the boundary fermions: 

q{n) = P+^{n)i + P^i){n)N, 

q{n) = ^{n)NP+ + ^(n)iP_. (II.3) 

The axial transformation is defined to rotate the right and left mode of the Dirac fermion 
oppositely. In the domain-wall fermion the massless mode is localized at the boundary of the 
flavor space with the right and left mode in the different side. The axial transformation of 
the domain-wall fermion is given by the flavor dependent vector like rotation of the fermion 
field which transform the different boundary fermions with opposite charge p|. 



5Xn), = e + - - zAXn),, (II.4) 
S^ii^n), = -e (^^ + 1 - ^^(n),A^ (II.5) 

where A" is a SU(A^/) generator of the rotation where Nj is a number of "real" quark flavors, 
and e(s) is a step function, 

^(■') = { -i: : < a ("■«' 

This transformation acts on the quark field as a usual axial rotation, 

S^qin) =iX'^^q{n), S^q^n) = iq{n)^^\\ {11.7) 

The corresponding axial current is given by the sum of the conserved vector current in the 
Wilson fermion over the flavor index, 

K(^) = Y.^{^ + \-s)j;{n,sl (II.8) 
J^M ^) = \ (^(^)« (1 + 7/.) U,{n)X-ij{n + /i), - i,{n + (1 - 7^) Ul{n)\^ij{n)\ . (II.9) 
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Unfortunately this transformation is not a symmetry of the action and therefore the axial 
current is not conserved. The divergence of the axial current becomes 

V^A'^in) = 2mg(n)A'^75g(^) - 2X'^(n), (11.10) 

where is a backward derivative and X" is the explicit breaking term characteristic for 
the domain-wall fermion, 

X"(n) = V'H|:P+AXn)^+i - ^{n) N^^P^X^^in) n. (11.11) 

In Ref. PI it is argued that the proper axial Ward-Takahashi identity is realized in — cx3 
limit if we consider the Green's functions with physical quark operators, 

{A;{n)0) - 2m {q{n)X'^^q{n)0) + i {6''0) = 0, (11.12) 

where operator O is made of the quark field ( [11. 3|) . The reason for this is that the explicit 
breaking term X'^ and the physical operator O is separated by N/2 in the flavor space and 
the Green's function (X"(n)(9) is suppressed by a factor e^"^/^, where a is some positive 
number. 

In this paper we set O = q{y)q{z) and calculate the one loop correction to the Ward- 
Takahashi identity. 



III. ONE LOOP CORRECTION TO WARD-TAKAHASHI IDENTITY 

In this section we consider the following Ward-Takahashi identity with O = q{y)q{z), 

= (A;(n)g(y)g(z)) - 2m ((g(n)A'^75g(ri)) q{y)q{z)) + 2 {X%n)q{y)q{z)) 

+Sn,y\^l5 {q{n)q{z)) + {q{y)q{n)) ^,\\ (111.1) 

Here the axial current for domain-wall quarks is given by a sum of the conserved vector 
current of the Wilson fermion system. One may wonder how this "vector" current turns 
out to be the continuum axial vector current '07^75''/' iii the a ^ limit. To see this we 
flrst study the continuum form of the Green's function {A'^^{n)q{y)q{z)^ at tree level and in 
N —>■ oo limit before the one loop calculation. In momentum space the tree level Green's 
function is given by Fig. [l| and has the following form in the continuum limit, 

{A^,{q)q{k)qip))^ = {qi~k)^ik, s)) 7^e(s, t) (^(-p, t)g(p)) , (111.2) 
where e{s,t) is a diagonal matrix deflned by 

<s,t)=e[^j + ^-sys, (II1.3) 

and {q{—k)il){k, s)^, {ip{—p, t)q{p)) are the external quark line propagator 0. A summation 
over the same flavor index is taken implicitly in this paper. Hereafter we will omit the 
SU(A^/) generator A'^ for simplicity. 
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We first notice the fact tliat tlie quark external line propagator is written as a product 
of the physical quark propagator and some damping factor in the continuum, 



Wn 



i|( + (1 — WQ)m 

1 — Wq 



2^ + (1 — WQ)m 



(III.4) 



(III.5) 



We have expanded the propagator to the next to leading order in the quark external momen- 
tum and mass, both of which give the leading order contribution in the one loop corrections. 
Now we can use the formula 

{w^'P^ + w^,-'P^) e(s, t) = {-w^-'P- + 1P+) = 75 {wfP- + ^0"'^) , (ni.6) 
{wt'P+ + w',-'P^) eis, t) = {-wt'P+ + w'o''P-) = -75 {wt'P+ + ^0"'^-) , (HI.T) 

which is valid at ^ 00 limit neglecting Wq^'^ terms. For this formula 75 factor appears 
and the leading order term of Green's function becomes 

AUMmP)), = 27.75 ., \, • (III.8) 

^ '0 + [1 — w^jm 1 — Wq ip + [1 — WQ)m 



This is a proper form of the Green's function in the continuum except for the overall factor 
1 - wl 

The tree level Green's function with the explicit breaking term vanishes in the continuum. 
{X^{q)q{k)q{p))o = E {^i-mk, s)) {p+6,.6,m^, - P-6,m^,6,.) 

s,t 

xm-p,mp))- (ni.9) 

It is easily shown that this is suppressed by Wq^'^ in — > 00 limit because of the damping 
factor in (|TL|) and (pL5| ). 



Now we calculate the one loop corrections to each terms in the Ward-Takahashi identity 
( [111.1| ). We set the external quark momentum and the quark mass to the physical scale, 
which is much smaller than the cut-off. We make an expansion in these variables and pick 
up the leading order terms which are relevant for renormalization. The next to leading terms 
are higher order errors in the lattice spacing. 

The one loop correction to the quark propagator and the pseudo scalar density vertex is 
already evaluated 0-0. The one loop level full quark propagator and the Green's function 
with pseudo scalar density become 



+ (1 - w^)Z^Z^^m 



wl) Z^Z2 



+ {l-wl)Z.^Z-^m^ ^ '''ii> + {l-wl)Z^Z-^m' 

(III.ll) 
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where ^2, Zm-, Zp^ and Z^ are lattice renormalization factors for the quark wave function, 
the quark mass, the pseudoscalar density operator, and the overall factor Wq, respectively. 
(A° is still omitted.) Their definitions are given in Ref. [Q. 

A one loop contribution to the Green's function (^A'^^{q)q{k)q{p)^ is given by the diagrams 
in Figs. ^ and ^ with vertices Ai ■ e which come from the axial vector current. The one loop 
diagrams in Fig. ^ contribute to the current vertex correction, 

/ , / \vertex 1 - "^n ^ ~ "^l 

{Al{q)q{k)q{p)) = —— \ ^ ^ ^^^^^ 

(III.12) 

where T^^^ is obtained by multiplying the vertex correction with the external propagator 
damping factor, 

.7.75 = Ur-'P. + wV'P-) [Vf + f « + f + m {s,t) 



Wo 



X 



(i.r-*)p„+u;r^)p+). (III.13) 



f '•'•''s are contributions from each diagrams represented by the massless fermion propagator 
Sp, fermion-gluon vertex Vi^, gluon propagator GJJ^ and axial current vertex A^^\l) given 
in the appendix. 





F 4\iT. ^i:(0> l)SF{l)s,uAf{l)e{u, u')Sp{l\,,Vl{-l, 0)G5(0, 


(III.14) 




r 7^ E ^1^(0, l)SF{l)s.uAf\l)Gt^{l)e{u, t), 


(III.15) 




£ E u)A(^>{l)SAlUVt{-l, 0)G^t(0, 


(III.16) 






(III.17) 



Here we have set the external quark momentum k^, and the quark mass m to be zero in the 
loop because the terms dependent on them are 0{a) and do not affect the renormalization. 

We can easily see that e(s,t) is multiplied directly by the damping factors in ( |111.13[ ) 
for r^^), r^f) and r|f^ and we can use the formula ( |111.7|) . For V^^^ we have a relation for 
massless fermion propagator. 



P+ iciw^-' + C2e-"(^-*)) + P_ (d^wl^ + rf2e-"(*-^))l , (III.18) 



+ 

where a and the coefficients of each term are explicitly given in the appendix. By neglecting 
^-aN/2 Qj^fjgj^ terms together with Wq^"^ we can make use of the formula similar to (pi.7| ). 
Finally the effective vertex becomes 
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(1) 



with 



- dH 

- dH 



up 



. (27r)4 t 

- dH 



- (27r)4 ^ 



-n (27r)^ 



V 

5:A«»(/)5^(/),,yit(-/,0)Gj(/), 



(III.19) 

(III.20) 
(III.21) 
(III.22) 
(III.23) 



Here we have used the fact that r^^*^ are proportional to 7^ and therefore anticommute with 

75- 

The half-circle and tadpole diagrams in Fig. ^ contribute to the wave function renormal- 
ization, 



iff + [1 — w^jm 



(N~s) - 



{wr^'^P. + wt'^P^) - {wt'^P. + 



Wq Jr-4 



Wo 



+ (1 — WQ)m 



xnKm)s,tSF{k)t,t'lAt'M {w^o'''^^P.+wi;'-'>P+) 



Wo 



1-w-^ 



tnwt'^P+ + w'o'-'"P 



l-w' 



1-wr 



2^ + (1 — Wo)m 



(N-s) 



+ (1 — 1^0)?^?' 



^111.24) 



We set the external quark momentum and mass to physical scale and extract the leading 
order terms in these variables. The fermion self-energy is the same as that is given in our 
previous paper . The difference is the existence of the fermion propagator in the internal 
line, which is expanded as follows together with the external line damping factor, 



{wr'^P^ + wt''>P^)SAp)s,t 



+ {1 — Wo)m 



w^-'P+ + wl'P^) + mwo {wl-^P+ + w^-'P^ 



{wl'P+ + 



w^'^P^ ) + mwo 



{<-'P. 



<-'P- 



1 



+ (1 — Wo)m 



(III.25) 
(III.26) 
(III.27) 



Here we notice that the leading and the next to leading order terms give the same order 
contribution. 
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With this expansion we have 



1 1 - w^, 

ii^ + (1 — WQ)rnP ^'^^ i]/) + (1 — tvDm 



Ey{p, m) 



l-wl 



+ (1 — Wo)m 
where E„ is defined by 



l-wl 



ilf, + {1 — ^ iji + (1 — Wo)rn' 



i|i + (1 — Wo)m 
1 



X 



The quark self-energy E^j is expanded in terms of and m, 



E(p, m)st = ^{0)st + ^Itt^^Pm + 



^ ^^^m + 0{p'^,m'^,pm), 



dipn dm 
and E^ becomes 



where 



(III.28) 



(111.29) 



^i^^ = K-^p+ + -r^p-) K-*p- + <^p+) , 

Summing up all the contribution we have one loop level full Green's function, 



(III.30) 



(III.31) 



(III.32) 

(III.33) 
(III.34) 



{A,{n)q{y)q{z))T' = {A,{n)q{y)q{z)), + {A,{n)q{y)q{z))-'^'' + (A,(n)g(i/)g(z))r™ 
{l-wl){Z.^Z2)-^ Za {I - tvDjZ^Z^Y^ 

+ wDZ^Z^'m 1 - wf^'^'t]/) + (1 - wDZ^Z^^m' 



(III.35) 



where 



Za^I + z!^^ +T^2\ 

Z2 = l + Zf\ 

z^ = i-zw, 

z-' = i-z^^^ + z^^\ 



(III.36) 
(III.37) 
(III.38) 
(III.39) 



In the following we will show the relation 
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9E(0) 



St 



and consequently Zg^^ = —T^\ We start by writing m) explicitly 

I 

up 



St 



n (27r)4 



+ 



A derivative with is given by 

as(o) 



v2Z{-p,p)G:im,t. 



where 



j(i) 



/-(2) 



/-(3) 



"^ '^ 



n (27r)4 



E 



up 



dip. 



(277)4 



Ev^i:(o,o^f(Om 



(27r)4 azp^ 

For the second equality in each equation we use the relation 

dSp{p + l),^t 



dipp 

dV,l{-p,l + p)s,t 



dipp 
V^p{-l-p,p)s,t 



dipp 
V2ti-P,P)s,t 



dipp 



Pm=0 



Pm=0 



Pm=0 



Pm=0 



Sf{1)suA^{0)Sp{IU, 

A';^%i)s,js,t, 



fj, J ^pu^s,t- 



The Eq. ( [1II.47| ) is given by differentiating the definition of the propagator 
"27m sin + W^ip)P+ + W-ip)p] Spip) = 1 



(III.40) 



(III.41) 



(III.42) 



^i(-^0)Gf,(/) = -4°)(.,t), (III.43) 



SFii)s,tVli-i,o)GZ{i) = -r«(^,^), (111.44) 



(III.45) 
(III.46) 



(III.47) 
(III.48) 
(III.49) 
(III.50) 

(111.51) 



and multiplying Sf{p) from the left. 

Equalities in ( |111.43|) to ( [111.46|) prove eq. ( |1II.40| ), which leads to = 1- This means 
that the axial vector current does not receive renormalization. 
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The one loop correction to the exphcit breaking term is given by the first diagram in 
Fig. ^ and those in Fig. § whose current vertex is now replaced by 



= \ P+6s,N/2^t,N/2+l - P~Ss,N/2+A,N/2) A". (III.52) 

In the diagram of Fig. ^the vertex X" is directly multiplied by the damping factor in external 
quark line and gives the suppression factor Wq^'^. A contributions form the remaining part 
of the diagram is 0{1) at most and the total contribution, which is a product of the two, 
vanish in X ^ cxd limit. For the first diagram in Fig. |^ the damping factor appears in 
the contracted fermion propagator, whose explicit form is given in the appendix. Since the 
vertex X"- is localized around s,t N/2, the damping factor produces suppressions such 
as w^^"^ or e~"^/^, and therefore the contribution from Fig. ^ again vanishes in the large N 
limit. 

Consequently the renormalized axial Ward-Takahashi identity is satisfied without fine 
tuning, 

n- f-h _L ^ I (1 - wl)Z^Z2 

~ ''^^'hli + (1 - wDZ^Z^^m^ ^^^^^'ri> + (1 - wl)Z.^Z;^^m 

~ \H + {l~ll)Z^Z-^m^^ - ^-^'"^^'WW^^^i^B^ 
^ {^-OZ^Z, {l-wl)Z^Z, ^ 

^ + (1 - OZ^Z-^m + + (1 - wDZ^Z-^m^'^ ' ^'''-"'^ 



together with conditions such that = 1 and Zp = Z, 



-1 

m 



IV. CONCLUSION 

In this paper we have calculated the one loop correction to the axial Ward-Takahashi 
identity and a renormalization factor for the "conserved" axial vector current in domain-wall 
QCD. Starting from the Green's function with the operators constructed from the physical 
quark field, we find that the axial Ward-Takahashi identity holds exactly without fine tuning 
at this order of the perturbation theory, and consequently the axial vector current defined 
by Furman and Shamir receives no renormalization, Zj^ = 1. As discussed in Ref. 0, one 
may expect that this property holds in all orders of the perturbation theory as long as zero 
modes exist. So one should use this almost conserved axial current to extract the pion decay 
constant in numerical simulations. 
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APPENDIX A. FEYNMAN RULES 



The fermion propagator is defined by 
SF{p)st={-i7^^smp^ + W-) ^GR{u,t)P+ + {~t^^smp^ + W+)^^GL{u,t)P^ (IV.54) 



with 



f-W 



W, 



s,t 



-W 



w. 



/-W 

1 



m 



1 

V m ~W J 

W = 1-M- r ^(1 - cosp^) 



s,t 



-W 



1 -wj 



(IV.55) 
(IV.56) 



GRis,t) = - 



(1 - m^) (l - VTe"") e"(-2^+^+*) - (1 - m^) (1 - VTe") e-"(^+*) 

+ v4e-"l^-*l, (IV.57) 



-2iy sinh(a)m (e°(-^+^-*) + e°(-^-^+*) 



t) = ^ [-(1 - m^) (1 - VTe") e°(-2^+-+*-2) _ (i _ m^) (i _ i^e~") e 



«{-s-t+2) 



cosh (a) 



-2W^sinh(a)m (e°(-^+^-*) + e"(-^-^+*) 

1 + W^' + E;.sin^p^ 
2iy 



+ Ae 



A 
F 



2Wsmh{a) 

1 - e^'W - (l - We"") . 
The physical quark propagator is given by 

—ij^ sinp^ + (1 ~ We~") m 



{q{-p)q{p)) 



- (1 - e'^W) + m2(l - We-'^) ' 
The gluon propagator can be written as 



(IV.58) 
(IV.59) 
(IV.60) 
(IV.61) 

(IV.62) 



Gtip) 



4sin^ p/2 



4 sinp^/2 sinpy/2 
4sin^p/2 



(IV.63) 



where sin^p/2 = sin^p^/2. The fermion-gluon interaction vertices which are relevant for 
the one loop calculation are given by 



Kt,{k,p)st = -igT'^i-ff.cos^i-k^+p^) - ir sinh-k^ + p^)}Sst, 



(IV.64) 



2^ M 2 
V2Uk,p)st = ^^?^^{T^T''}{^7Msini(-fc^ - r cos^i-k, + p,)}5,J,,. (1V.65) 

The current vertices used in Eqs. ( |111.14| ) to ( [111.17| ) are given by expanding the axial 
vector current in terms of the gauge field. Only three kinds of them are relevant for one 
loop calculation, 
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Ho = liJ. COS - ir sin Z^, 



fl' ( 7m sin ^ +ircos^ J T"", 



(IV.66) 

(IV.67) 



A{2)ab 



9^ 



arpb 



(IV.68) 

The exphcit form of the contracted propagator with the damping factor is given by 



where 



( 



Sf{p) 



+ {{wo - W{p)) w'o~^Gr{s, t) + mGR{N, t)) P+ 
+ {(wo - W(p)) w^-'Gl(s, t) + mGL(l, t)) P_, 

SF{p)t,s {;wi''-'^P^ + wt'^P.) 

= {P+wt'^GLis,t) + P,wil'-'^Gnis,t)) i-t^,smp,) 
+P+ {wi''-'^GL{s, t) {wo - W{p)) + mGLit, 1)) 
+P_ {wf^Gnis, t) (wo - Wip)) + mGnit, N)) , 



w',-'Gr(s, t)^j [-{1 - m') (1 - We'') e 
-2iy sinh(Q;)me"" 



-2a 



- Woe" 

-a{N-t) 



,-a(t-l) 



woe- 



l^^e-°(*-'-Hl-»or--')-2»o 



(f-i) 



a\] sink a 



1 + Wq — 2wo cosh a 



A 



w^-^Gn{s, t) = - [-(1 - m') (l - We 



1 



1 — WqC' 



-a{N-t) 



-2W^sinh(Q;)me-"Y 
e-"(^-*)(l - Woe 



-a(t-l) 



3— a 



+A 

A 



Woe 

^"'^ — 2wo^ ^^wq sinh a 



1 + Wq — 2wo cosh a 



w',~'Gl{s, t) = - -(1 - m^) 1 - We-" 



1 



— 21^sinh(Q;)me" 



1 



1 — woe~ 



1 — Woe " 

_ a(Ar-t) 



wt'GL{s,t)^^- 



+A 

A 



-a(t- 



"^^(1 — Woe ") — 2^0 ^Wosinh 
1 + Wo — 2wo cosh a 



a 



-(1 - m') (1 - We") 
1 



-2W sinh(a)me" 



1 — Woe " 



1 - Woe-°= 



-a{N-t) 



(IV.69) 



(IV.70) 



(IV.71) 



(IV.72) 



(IV.73) 
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A 



GR{N,t) = - 



1-\-Wq — 2wq cosh a 
(1 - m') (l - We"') e-"(^-*) 
-2PFsinh(a)me-"e~"(*-^)] + ^e'^^^"*), 



— 2Vr sinh(Q;)me 



a-a{N-t) 



^^g-a(t-l)_ 



(IV.74) 



(IV.75) 



(IV.76) 
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FIG. 3. One loop diagrams which contribute to the quark wave function. 
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